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Abstract. 

We present new observations of the center-to-limb variation of spectral lines in the quiet Sun. Our long-slit 
spectra are corrected for scattered light, which amounts to 4-8 % of the continuum intensity, by comparison with 
a Fourier transform spectrum of the disk center. Different spectral lines exhibit different behaviors, depending on 
their sensitivity to the physical conditions in the photosphere and the range of depths they probe as a function 
of the observing angle, providing a rich database to test models of the solar photosphere and line formation. We 
examine the effect of inelastic collisions with neutral hydrogen in NLTE line formation calculations of the oxygen 
infrared triplet, and the Na I A6160.8 line. Adopting a classical one-dimensional theoretical model atmosphere, 
we find that the sodium transition, formed in higher layers, is much more effectively thermalized by hydrogen 
collisions than the high-excitation oxygen lines. This result appears as a simple consequence of the decrease of the 
ratio Nh/N c with depth in the solar photosphere. The center-to-limb variation of the selected lines is studied both 
under LTE and NLTE conditions. In the NLTE analysis, inelastic collisions with hydrogen atoms are considered 
with a simple approximation or neglected, in an attempt to test the validity of such approximation. For the 
sodium line studied, the best agreement between theory and observation happens when NLTE is considered and 
inelastic collisions with hydrogen are neglected in the rate equations. The analysis of the oxygen triplet benefits 
from a very detailed calculation using an LTE three-dimensional model atmosphere and NLTE line formation. 
The x 2 statistics favors including hydrogen collisions with the approximation adopted, but the oxygen abundance 
derived in that case is significantly higher than the value derived from OH infrared transitions. 
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1. Introduction 

Abundance analyses of late-type stars commonly rely on 
the assumption of Local Thermodynamical Equilibrium 
(LTE) and a plane-parallel homogeneous structure. Even 
though in many cases the derived abundances may be ac- 
tually similar to the real photosphcric abundances, there 
are well-known instances when departures from LTE and 
inhomogcneities can introduce significant systematic er- 
rors. Consequently, efforts have been directed toward im- 
proving the modeling techniques. One of the most signifi- 
cant obstacles in the way of performing reliable Non-LTE 
(NLTE) calculations is the incompleteness of the necessary 
atomic data. Because the calculations tend to be involved, 
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and rely on approximations to account for rates that have 
never been measured in a laboratory or cannot be calcu- 
lated reliably, it is often the case that results cannot be 
directly accepted without extensive testing and thorough 
comparison with observations. 

Common tests that have been put into practice are 
spectroscopic observations of profiles of lines with well- 
known damping constants, abundance determinations for 
stars in a cluster (expecting all members to show the same 
chemical composition) or a solar abundance analysis for 
non-volatile elements (anticipating that the photospheric 
abundances and those found in chondrites are identical). 
The spatially resolved solar disk offers the possibility to 
survey a range of formation depths for the continuum, 
and also for any given spectral line. The center-to-limb 
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variation of spectral lines constitutes a particularly in- 
teresting test for elements like oxygen, which is depleted 
in meteorites (see e.g. Sedlmayr 1974; Kiselman 1991). 
Furthermore, this type of observations may turn very use- 
ful to distinguish between the thermalizing effect of inelas- 
tic collisions with electrons and hydrogen atoms, as the ra- 
tio N e /NH increases by more than one order of magnitude 
between logrsooo = ~ 2 and 0. 

Inspection of the available literature reveals a fla- 
grant scarcity of high quality observations of the center-to- 
limb variation of line profiles. Early work did not benefit 
from digital detectors (e.g. Miiller & Mutschlecner 1964; 
Miiller, Baschek, & Holweger 1968). The extended study 
by Balthasar (1988) includes more than one hundred lines, 
but only reports equivalent widths and a few parameters 
related to the amount and shape of the line asymmetry. 
Equivalent widths, however, neglect much of the available 
information in a line profile. Equivalent widths ignore, for 
example, whether changes in the damping wings or the 
core of a line dominate the line strength variation, or even 
if they cancel each other to produce an equivalent width 
nearly independent of the position on the disk. More re- 
cent work exists, but it is limited to a small number of 
spectral regions (e.g. Ambruoso et al. 1992; Brandt & 
Steinegger 1998; Grigoryeva & Turova 1998; Langhans & 
Schmidt 2002; Stenflo et al. 1997). 

We have obtained observations of a number of key 
lines for which reliable damping constants are available. 
In Section [5] we report our measurements and the data 
reduction, in particular a procedure to subtract the scat- 
tered light. In Section|21we study the cases of the infrared 
oxygen triplet and the Na I A6161 lines in the context 
of a classical model atmosphere. In Section 0J the case 
of the oxygen triplet is reanalyzed with a time-dependent 
three-dimensional model of the solar photosphere. The pa- 
per concludes with some reflections about the results and 
suggestions for future work. 

2. Observations and data reduction 

Solar observations of the center-to-limb variation of sev- 
eral spectral lines were carried out in October 22-23, 1997, 
with the Gregory Coude Telescope (GCT) and its Czerny- 
Turner echelle spectrograph (Kneer et al. 1987; Kneer 
& Wiehr 1989) at the Observatorio del Teide (Tenerife, 
Spain). This telescope was moved to Tenerife and refur- 
bished in 1985, after more than 20 years of operations in 
Locarno (Switzerland). It was dismantled in 2002 to leave 
room for new instrumentation. The spectrograph was op- 
erated in orders 6-9 with a slit width between 50 and 150 
/xm (0.4-1.3 arcsec) achieving an estimated FWHM (= 5X) 
resolving power R = X/SX in the range 57,000-240,000. 

We secured spectra for 8 spectral setups in 6 different 
positions across the solar disk, as summarized in Table 1. 
A ninth setup was centered at about 7610 A, to obtain an 
estimate of the amount of scattered light from saturated 
telluric O2 lines. The exact slit location was always chosen 
to avoid active regions. Ten consecutive equal-length ex- 



Table 1. Observations. 



Central A R' S (Set. light) max(6>) 

(A) < a > (deg) 



5245 


56000 


0.054 


75 


5300 


176000 


0.046 


80 


6122 


77000 


0.065 


80 


6162 


148000 


0.066 


75 


6166 


77000 


0.066 


75 


6200 


206000 


0.080 


75 


7585 


176000 


0.048 


80 


7770 


86000 


0.058 


80 



K 1 is the resolving power measured relative to that of the 
FTS spectra; see text 
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Fig. 1. Correspondence between the location and extent 
of the slit for each position in /i = cos 9 and the angle 
on the sky along the observations line from the center of 
the disk lo = sin where fl is the solar angular radius 
which is assumed to be 959.5 arcsec. The positions of the 
center of the slit are marked with asterisks. For position 
# 6, part of the slit was outside the solar disk; the center 
of the illuminated part of the slit is located at the tip of 
the arrow, and the right limit is marked with a vertical 
line at lo = 17. 

posures were obtained at each position for each setup. The 
exposure time varied, depending on the setup, between 0.5 
and 2.0 seconds. Positions #1 to #5 were always at helio- 
centric angles 9 = 0, 15, 30, 45, and 60 degrees (^ = cos 9 
= 1.00, 0.97, 0.87, 0.71, and 0.50) along a straight line 
crossing the center of the solar disk. Position #6 was also 
selected along the same direction, sometimes at 9 =75 de- 
grees and others at 80 degrees = 0.26 or 0.17). The 
slit length covered approximately 205 arcsec, but the field 
was truncated by the CCD size to w 160 arcsec. In the 
last position, part of the slit was outside the solar disk, 
and therefore when the center of the slit was at 9 =75 and 
80 deg, the center of the illuminated slit was at 70 and 72 
deg, respectively, or approximately fi ~ 0.32. Fig^shows 
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Fig. 2. Contour plot showing the reduced x 2 , normalized 
to unity at its minimum, as a function of the relative re- 
solving power and the amount of scattered light for the 
fi = 1 setup centered around the Ca I A6162 line. 

the slit extent both in angular size on the sky and in fi 
for each position, assuming a solar angular radius of Q = 
959.5 arcsec (Allende Prieto et al. 2003a). In what follows 
we neglect limb-darkening within the slit length, assigning 
any observed intensity to the location of the center of the 
illuminated part of the slit. 

No calibration lamps were available. The detector 
was a 1024 2 CCD with negligible dark current (1.4 
ADU/sec/pixel). The CCD frames were bias subtracted, 
and rotated. The rotation corrects small tilts between the 
direction of the spectral dispersion and the CCD that 
ranged between 0.5 and 0.7 degrees, as derived from the 
drift of the centers of several spectral lines in the spa- 
tial direction. Pixel-to-pixel variations and weak fringes 
were automatically smeared out in the spatially-averaged 
spectra at each position. Wavelength calibration of the 
disk-center spectra was carried out by fitting a 1st to 3rd 
order polynomial (depending on the number of available 
spectral lines) to the wavelengths of the line centers, as 
measured in the Brault & Neckel atlas (1987; see Neckel 
1994). The wavelengths were adopted from Allende Prieto 
& Garcia Lopez (1998), or measured afresh when missing 
from their list. The same dispersion solution was applied 
to the observations at other positions. Finally, a velocity 
shift relative to the center-of-the-disk spectra (accounting 
mainly for solar rotation, Earth's motion, and the limb 
effect) were determined by cross-correlation and applied 
to the spectra at other positions. Although all the obser- 
vations for a given spectral setup were acquired consecu- 
tively, the time interval between the observations at the 
disk center and at the limb was generally significant, and 
therefore we deem our wavelength accuracy insufficient to 
study the limb effect. 

The spectrum of O2 lines at about 7600 A (A band) 
confirmed the suspicion that scattered light was substan- 
tial. The presence of scattered light is also apparent when 
comparing our spectra with the Fourier transform (FT) 
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Fig. 3. Comparison of the spectra at the center of the 
disk in the region around Ca I A6162.2 from Brault & 
Neckel (1987; convolved with a Gaussian; solid line) and 
the Gregory Coude Telescope (after removal of scattered 
light; filled circles). 

spectrum at the center of the disk in the Brault & Neckel 
atlas. Fourier transform spectrographs (FTS's) are not 
susceptible to scattered light in the same sense as grat- 
ing spectrographs, but non-linearities in the detector can 
cause systematic errors in the zero-point of the intensity 
scale (see the discussion in Kurucz et al. 1984). Scattered 
light leaves a characteristic distortion on the spectrum. On 
visual inspection, for a given spectral resolution, the cores 
of the lines appear filled. This effect is noticeably different 
from a reduction in spectral resolution. As the spectrum 
of the quiet Sun averaged over a large area is believed to 
be extremely stable, two instrumental factors are mainly 
responsible for the differences between the FT spectrum of 
Brault & Neckel and our center-of-the-disk spectra: resolv- 
ing power, and the amount of scattered light. Therefore, 
based on the well-supported expectation that the amount 
of scattered light in the FTS atlas is much less that in 
our GCT spectra, we used the FTS data as a template to 
correct the scattered light 1 . 

For each of the setups, we compared our center-of-the- 
disk spectra with the FTS atlas to derive the amount of 
scattered light and the resolving power (R' relative to the 
FTS atlas 2 ) that would lead to the best agreement be- 
tween the two. Our analysis was based on modeling the 
instrumental profile of the GCT spectra as a Gaussian, 
and the scattered light as a constant fraction of the con- 
tinuum flux for each setup. The comparison required the 
GCT spectra to be first continuum corrected. This was 
accomplished by using a 6th-order polynomial and a se- 

1 The O2 lines of the A-band, which are presumed opaque, 
show ~ 0.5 % of the pseudo-continuum intensity in the FT 
spectrum, while our spectra indicate > 2 %. 

2 The relative resolving power R' can be approximately con- 
verted to absolute resolving power knowing that Rfts — 
4 x 10 5 , and R = (1/R' 2 - l/i?| TS ) " 1/2 . 
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Fig. 4. Center (filled, circles) and limb (solid) spectra in 
the region around Ca I A6162.2 (top panel) and around 
Co I A5301.0 (bottom panel). 



ries of clipping iterations. Then a search was performed 
to determine the best-fitting values for the FWHM of the 
Gaussian representing the instrumental profile, and S, the 
amount of scattered light expressed as a fraction of the 
pseudo-continuum intensity. 

To avoid a bias due to the original normalization, the 
procedure was iterated until SR'/R' < 10~ 3 . Only 3-4 
iterations were required to comply with the convergence 
condition. Table 1 provides the results. Figure [21 shows 
the contours of the reduced x 2 (normalized to 1 at the 
minimum) for the setup centered around the Ca I A6162 
line. The signal in the GCT spectra has the potential to 
exceed a S/N of at least 1000. Systematic uncertainties are 
more difficult to assess. Systematic errors are expected due 
to shortcomings of modeling the GCT instrumental profile 
as a Gaussian, and the scattered light as independent of 
wavelength (for a given setup) , perhaps with an additional 
contribution from real variations of the photospheric solar 
spectrum with time. 

We estimate that our derived values of the resolving 
power are good to about a few percent, and that the scat- 
tered light is constrained to within ~ 1 %. Both the rel- 



ative scattered light, and the polynomial representing the 
pseudo-continuum were applied to correct the spectra at 
the center of the disk (from which the values were deter- 
mined) and also to the spectra at other positions. Fig. 
13 shows the FTS spectrum smoothed to the resolution of 
the GCT data, and the GCT spectrum after correcting the 
scattered light. We empirically find that the observed pro- 
files are accurate to ~ 1 %. This figure is also supported by 
an rms difference of ~ 0.7 % between two corrected spec- 
tra of the 6122 A region at the center of the disk obtained 
on the two consecutive observing dates. 

Fig. 01 shows the change from the center to the limb of 
the spectral region discussed in the previous figures (up- 
per panel), and for the setup centered at w 5300 A (lower 
panel). Many differences are noticeable. Some lines only 
undergo subtle changes, often becoming slightly broader 
towards the limb, but other features' wings turn less pro- 
nounced. Some lines become both wider and stronger to- 
wards the limb (e.g. Co I A5301.0). The variations in be- 
havior among different lines reflect their different sensitiv- 
ities to the physical parameters controlling the line forma- 
tion and, as earlier recognized in the 60's, contain precious 
information on the atmospheric structure. The GCT spec- 
tra have been made publicly available through the CDS. 
In the next section we will explore the possibility of quan- 
tifying the thermalizing effect of inelastic collisions with 
electrons and hydrogen atoms in the populations of neu- 
tral oxygen and sodium through the center-to-limb behav- 
ior of spectral lines. 

3. Center-to-limb variation of spectral lines and 
departures from LTE 

Line profiles sample a range of atmospheric depths that 
changes depending on the position on the disk. An analy- 
sis of the resulting changes in the line shapes is, to a first 
approximation, independent of the /—value of the tran- 
sition, the chemical abundances, and the damping con- 
stants. A zeroth-order analysis can be based on simplified 
LTE calculations. Such modeling already reproduces the 
basic behavior for many spectral lines; their strengthening 
or weakening from center to limb due to the changes in the 
equilibrium populations with atmospheric depth. A closer 
look at the observed line profiles will reveal the intrica- 
cies of the line formation, such as departures from LTE, 
the presence of inhomogeneities, or the lack of complete 
frequency redistribution. 

We explore first the possibility of discerning the ther- 
malizing effect of collisions with electrons from collisions 
with hydrogen atoms by means of an analysis based on 
one-dimensional model atmospheres and NLTE line for- 
mation. When dealing with high S/N and resolving power 
spectroscopic observations, an analysis based on static 
model atmospheres yields a very poor match between cal- 
culated line profiles and observations. This is mainly the 
result of neglecting surface inhomogeneities - granulation- 
which cause the observed lines to appear asymmetric. 
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Fig. 5. Response functions to the temperature at /i = 0.32 
and 1.00 for the core of the lines of the O I infrared triplet 
(solid) and the Na I A6160.8 (dashed) lines. 



Admittedly a handicapped approach, such a detour is con- 
venient and oftentimes didactic. 

We examine two interesting cases which are part of our 
observations: the oxygen infrared triplet, and a Na I line 
at 6160.8 A. Fig. shows the LTE response to temper- 
ature perturbations of the core of these lines. These re- 
sponse functions (see Ruiz Cobo &: del Toro Iniesta 1992, 
1994) are calculated for the continuum-corrected line pro- 
files, and therefore include the sensitivity of the contin- 
uum, which strengthens the lines when the temperature 
increases in the deepest layers. In this plot, a positive value 
represents an increase in the relative flux, and therefore 
a weakening of the line. When changing from the center 
to the limb, the O I lines map fairly well the layers in the 
range —0.5 < logr 50 oo < +0.5. The temperature sensitiv- 
ity of the core of the Na I line shows opposite reactions to 
changes in the temperature at different heights but this 
line is sensitive to higher layers than the triplet. 



(Kurucz 1993). The model atom and the calculations fol- 
low those in Allende Prieto et al. (2003a, 2003b). We 
have introduced the effect of collisional excitation and ion- 
ization due to neutral hydrogen with the prescription of 
Steenbock & Holweger, multiplying the rates by a cor- 
rection factor 5h- To make the test independent of the 
oxygen abundance, we have adjusted it for each of the 
following three cases: LTE, Sh = 0, and Sh = 1. A micro- 
turbulence of 0.9 km s _1 was adopted - a typical value 
for the disk-center (e.g. Blackwell, Lynas-Gray, & Smith 
1995). A Gaussian macro-turbulence of 2.0 km s _1 was 
used, as well a Gaussian instrumental profile consistent 
with the analysis in Section [21 (see Table 1), for a total 
FWHM broadening of 0.151 A. Collisional damping in- 
duced by neutral hydrogen was accounted for using the 
cross-section at 10,000 K by Barklem, Piskunov & O'Mara 
(2000), assuming a temperature dependence T 2 / 5 . 

Figure compares observed and calculated profiles 
for the different positions and the three cases considered: 
LTE, Sh = Q, and Sn = 1. The first three panels illustrate 
the agreement between observed and computed profiles at 
the extreme positions: \i = 1.0 (center) and 0.32 (near the 
limb). The fourth panel quantifies the relative agreement 
through x' 2 = XX .fobs - jsynth) 2 , which has been normal- 
ized to 1 at the minimum. To avoid obvious blending fea- 
tures, only the fluxes which are more than 3 % lower than 
the continuum level are considered in the evaluation of 
X 12 - As collisions drive the level populations toward LTE, 
the NLTE line profiles with Sh = 1 are an intermediate 
case between LTE and the Sh = profiles. As we antici- 
pated, the match of the observed profiles with a static one 
dimensional model is poor. The comparison suggests that 
the NLTE calculations perform much better than those 
assuming LTE, in agreement with previous investigations. 
We can also conclude that the effect of hydrogen collisions 
is quite limited in these lines (see also Nissen et al. 2002), 
but the best agreement is found when they are included, 
albeit the difference is marginal. 



3.1. The oxygen infrared triplet 

It has been recognized that departures from LTE are sig- 
nificant for these lines. The main effect is due to an in- 
frared radiation field weaker than Planckian in the line 
formation region (Eriksson & Toft 1979). Fortunately, the 
level populations of O I are relatively insensitive to the 
UV radiation field in solar type stars, which simplifies the 
modeling. Due to cosmological and galactic implications, 
oxygen abundances in metal-poor stars have been exten- 
sively discussed. When modeling the oxygen triplet lines 
in NLTE some authors opt to neglect collisions with neu- 
tral hydrogen (e.g. Nissen et al. 2002), while others prefer 
to adopt a recipe first suggested by Steenbock & Holweger 
(1984; see, e.g., Takeda 2003), based on Drawin's formula 
(Drawin 1968), scaled by an empirical factor(s). 

We have evaluated the NLTE populations for O I 
with a solar model from Kurucz's non-overshooting grid 



3.2. Na I A6160.8 

Departures from LTE in Na I lines have also been exten- 
sively studied (e.g. Athay & Canfield 1969; Baumullcr, 
Butler, & Gehren 1998). The level populations are again 
highly independent of the UV radiation field. We have 
used the new model atoms and calculations we previously 
referred to for the oxygen case. A micro-turbulence of 0.9 
km s -1 was adopted, as in Section 13.11 but a lower value 
of the macro-turbulence than for the oxygen triplet was 
necessary to fit the profiles at the center of the disk, 1.3 
km s _1 , and therefore the total FWHM for the Gaussian 
broadening was 0.074 A (note the difference in resolving 
power, as summarized in Table 1). The effective princi- 
pal quantum number (n*) of the 2 S upper state is larger 
than 3 and therefore exceeds the range of the tables pub- 
lished by Anstee and O'Mara (1995). Barklem et al. (2000) 
derived damping cross-sections for other Na I transitions 
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Fig. 6. Center (at rest with the continuum normalized to 1) and limb (shifted to the red with the continuum normalized 
to 1 and subtracted 0.05) spectra for the O I infrared triplet lines. Observations (filled circles) and calculations are 
compared for three different models: LTE, NLTE/Sh = 0, and NLTE/Sh = 1. The bottom-right panel shows the 
variation of the quality of the fit as a function of the position of the disk: LTE (solid), NLTE/Sh = (dashed) and 
NLTE/Sh = 1 (dot-dashed). The x' 2 values have been normalized to 1 at the minimum. 



with similar energies sharing the same lower state 2 Py 2 - 
For those transitions the FWHM of the Lorentzian (colli- 
sional) component of the line absorption profile (per unit 
perturber number density at 10,000 K) log 7 was found 
to be in the range —6.86 to —7.23. Comparison with the 
observed profiles at the disk center showed that, indepen- 
dently of the chosen abundance, the best agreement is 
found for a value close to log 7 = —6.86, which was there- 
fore adopted. Accepting n* = 3 for the upper s state gives 
log 7 = —6.96. 

Fig. shows the profiles at the extreme positions 
(fi = 0.32 and 1) and the comparison with the synthetic 
spectra for the three cases considered. Again, we adjusted 
the Na abundance for each case in order to maximize the 
agreement between calculated and observed profiles at the 
center of the disk. It is evident that the sensitivity of the 
level populations involved in this Na I transition to colli- 
sions with H is higher than for the states connected by the 

I infrared triplet. This is most likely the result of the 
line formation and the departures from LTE taking place 
in this case in higher layers. The LTE response function 
to the temperature for the core of the O I triplet lines 
peaks at logrsooo — 0, while that for the core of the Na 

1 A6160.8 line does so at logT 50 oo — — !• Using the Van 



Regemorter (1962) and Drawin approximations, for a typ- 
ical transition of a few electron volts the ratio of the col- 
lisional rate coefficients qH/q c cx (U + 2)/U 2 e~ u J E±(U), 
where U = J5y/ (kT), and E± is the first-order exponential 
integral. Thus, qn/qe is a fairly flat function of the tem- 
perature or the optical depth in the solar photosphere. 
Between logr 50 oo — and —1, Nh/N 6 increases by an or- 
der of magnitude, making the effect of collisions with hy- 
drogen more noticeable in higher layers. Similarly to the 
case of the oxygen triplet, we find the best agreement with 
observations when departures from LTE are considered, 
but now the calculations that neglect collisions with hy- 
drogen grade best, and the differences between the Sh = 
and 1 cases are now more significant. Not even the NLTE 
Sh = model can match the wider profile of the sodium 
line that is observed near the solar limb. 

4. The O I infrared triplet lines in a 
time-dependent 3D model atmosphere 

Modern radiation hydrodynamics simulations offer a much 
more realistic representation of the lower solar atmosphere 
than classical static model atmospheres (Stein & Nordlund 
1998; Asplund et al. 2000). In the context of these, more 
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Fig. 7. Center (at rest with the continuum normalized to 1) and limb (shifted to the red with the continuum normalized 
to 1 and subtracted 0.05) spectra for Na I A6160.8. Observations (filled circles) and calculations are compared for three 
different models: LTE, NLTE/Sh = 0, and NLTE/Sh = 1. The bottom-right panel shows the variation of the quality 
of the fit as a function of the position of the disk as in Fig. |S| 



detailed, model atmospheres, surface convection and line 
broadening due to small- and large-scale velocity fields are 
naturally included in the calculations. This makes the use 
of ad hoc parameters such as macro- and micro-turbulence 
unnecessary, and at the same time improves substantially 
the agreement between observed and calculated spectral 
line profiles. 

The consideration of several more dimensions makes 
calculations, however, much more involved, which, with 
some notable exceptions, has largely limited NLTE stud- 
ies to ID modeling. However, O I is among the few ions for 
which the line formation in a 3D hydrodynamical model 
has been studied with a complex multi-level model atom. 
Asplund et al. (2004) have recently shown that it is nec- 
essary to consider both, departures from LTE and a 3D 
geometry, to find good agreement among different spectro- 
scopic indicators of the solar photospheric oxygen abun- 
dance, finally settling a long-standing problem. Here we 
will use similar calculations to revisit the center-to-limb 
variation of these lines, which were analyzed with a ID 
model in Section l3~T1 

The 3D solar simulation has been described in detail 
in previous papers, and the NLTE line formation calcu- 
lations are identical to those reported by Asplund et al. 
(2004). We proceed similarly to the ID analysis. Line pro- 



files for the O I triplet were computed for different incli- 
nations from the vertical axis, and averaged over azimuth, 
time, and horizontal position. The computed NLTE pro- 
files are the product of an average LTE profile derived 
from 100 snapshots covering 50 minutes of solar time, and 
the ratio of the NLTE and LTE line profiles derived from 
two snapshots. This procedure reduces the computational 
burden, while introducing negligible errors. The statistical 
equilibrium equations were solved for the two considered 
snapshots for two different cases: Sh = 0, and Sh = 1. 
Three different abundances were used in the calculations 
(loge(O) = 8.50, 8.70, and 8.90 3 ) and later the profiles 
were linearly interpolated to find the best match at the 
disk center in each of the three cases: LTE, Sh = 0, and 
Sh = 1. The calculated profiles were smoothed by convolu- 
tion with a Gaussian with a FWHM of 0.088 A to account 
for the finite spectral resolution of the observations. 

Fig. |H1 is the counterpart of Fig. for the 3D case. 
Admittedly, the computed line profiles are not perfect. In 
particular, significant discrepancies are noticeable around 
the cores of the lines. The reader should keep in mind 
that the lower- right panels of Figs. El and |H1 cannot be di- 
rectly compared, as the x' 2 is normalized to unity at the 

3 e(O) = 10 12 N (0)/N(H), where N represents number den- 
sity. 
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Fig. 8. Similar to Fig. for the 3D time-dependent simulation of the solar surface. 



minimum independently for each case: ID and 3D. The 
adopted abundances were chosen to optimize the overall 
fit of the /J, = 1 observations, and the agreement is con- 
siderably better in 3D than in ID near the line wings - 
even though macro-turbulence is included in the ID case 
and not in 3D. The refined 3D analysis strengthens what 
the ID study suggested. The calculated profiles are now 
close enough to the observations to justify an attempt to 
quantify statistically the agreement. Making use of the ac- 
curacy estimate empirically derived in Section[21 we adopt 
a ~ 0.01, and evaluate x 2 = 1/c 2 Ef = i (jobs - /synth) 2 for 
/U = 0.32 in the three cases. Using only the fluxes below 
3% of the continuum level to avoid the interference of ob- 
vious blending features between the lines, n= 135, and we 
find x 2 = 725, 195, and 92 and therefore significance levels 
of (< 10~ 81 ), 10~ 3 , and 1, for LTE, S K = 0, and Sh = 1, 
respectively. 



In the ID analyses we have avoided comparing abso- 
lute abundances, given that the effect of atmospheric in- 
homogeneities was neglected. Now that the triplet lines 
are modeled in more detail, it is tempting to examine the 
absolute abundances as an additional piece of evidence. 
The abundances found in each case to match the profiles 
at the center of the disk are loge(O) = 8.88, 8.66 and 8.72 
for LTE, 5 H = 0, and 5 H = 1, respectively. The 5 H = 
case provides the best agreement with the average value 



derived by Asplund et al. (2004) 4 : 8.66 ± 0.05 dex from 
forbidden and permitted lines of O I and from infrared OH 
lines. Thus, absolute abundance and % 2 statistics lead to 
apparently opposite conclusions in the case of the oxygen 
triplet. We should note that the average value derived by 
Asplund et al. (2004) includes permitted lines analyzed 
in a similar manner as in this paper with NLTE line for- 
mation and Sh = 0. A higher oxygen abundance from 
allowed lines is expected for Sh = 1, and in that case, for- 
bidden and permitted lines would suggest a higher abun- 
dance (8.67-8.70) than infrared OH lines (8.61-8.65) dex. 

We should stress that the necessary atomic data for 
the NLTE calculations were independently compiled for 
the ID and 3D cases, yet the qualitative agreement is ex- 
cellent. Quantitatively, the abundance corrections neces- 
sary to make the LTE and NLTE equivalent widths agree 
are also remarkably similar in ID and 3D. This was also 
found by Asplund et al. (2004) comparing the 3D calcu- 
lations with ID abundances based on MARCS model at- 
mospheres. The impact of departures from LTE and that 
of surface inhomogeneities on the abundances cannot be 
generally decoupled, but the situation for the OI triplet, 
where this is actually a good approximation, makes the 



4 Our derived oxygen abundance for the 5*h = case (8.66 
dex) is consistent with the abundances derived by Asplund et 
al. (2004) by fitting the profiles of the same lines in the solar 
flux atlas of Kurucz et al. (1984): 8.64-8.66 dex 
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conclusions from our ID analysis to be similar as those 
from a 3D study. 

5. Summary and conclusions 

The type of observations presented here, as exemplified 
in Sections [3] and provide a stiff test of theoretical cal- 
culations of atmospheric structure and line formation in 
a solar-type photosphere. The potential of center-to-limb 
observations to guide modeling was recognized early, but 
consistent fj, ^ 1 observations covering extensive spectral 
regions are still missing. We have explored how these ob- 
servations can constrain the effect of collisions with hydro- 
gen atoms on the rate equations for two particular cases. 

Our results for Na I A6160.8 fall in line with the con- 
clusion drawn from theoretical calculations by Barklem, 
Belyaev, & Asplund (2003) for excitation of another al- 
kali, lithium, by inelastic collisions with H atoms (see also 
the discussion by Lambert 1993). Limited by the accu- 
racy of the observations, and approximations involved in 
computing line profiles, we have only tested here the ex- 
treme cases of 5h = (no collisions with H), or Sh = 1 
(Drawin-like formula). Comparison with detailed NLTE 
calculations based on multidimensional time-dependent 
model atmospheres and improved solar observations free 
from scattered light may be able to answer whether or 
not the use of a simplified formula scaled by a constant 
factor (<< 1 for Na I or Li 1) is a useful approach. Our 
results and a quick inspection of the literature lead to the 
naive conclusion that such a factor would need to vary for 
different species. For example, the role of hydrogen colli- 
sions compared to those with electrons increases for atmo- 
spheres more metal-poor (or cooler) than solar, and Korn, 
Shi, & Gehren (2003) have recently found that Sh = 3 is 
needed to satisfy the iron ionization balance in several 
metal-poor stars. 

In the case of the oxygen infrared triplet, the solar 
observations are best reproduced considering inelastic hy- 
drogen collisions in the rate equations. However, the oxy- 
gen abundance obtained when hydrogen collisions are ne- 
glected (8.66 dex) is in better agreement with the abun- 
dance inferred from OH infrared lines (8.61-8.65 dex; 
Asplund et al. 2004) than the abundance derived from 
the NLTE 5 H = 1 analysis (8.72 dex). The higher S H = 1 
abundance is nevertheless still consistent with the values 
from forbidden lines (8.67-8.69 dex). The weighted aver- 
age of the abundances from permitted lines derived by 
Asplund et al. (2004) changes from 8.64 ± 0.02 dex to 
8.70 ±0.04 dex when inelastic collisions with hydrogen are 
considered with the Drawin-like formula. Therefore, this 
apparent inconsistency could signal that the abundance 
derived from infrared OH lines is more uncertain than the 
values obtained from a detailed analysis of forbidden and 
permitted O I lines. 

We should bear in mind that the effect of collisions 
with hydrogen on the statistical equilibrium is not the only 
important uncertainty for NLTE calculations in late-type 
stellar atmospheres. Electron collisions, for example, are 



only approximately included in our models, mainly due to 
the lack of reliable data for particular states. Our results 
should be confirmed by taking into account refined data as 
they become available (see, e.g., Zatsarinny & Tayal 2003). 
In addition, the electron density in high atmospheric lay- 
ers is highly uncertain even in the 3D hydrodynamical 
simulations, due to the assumption of LTE. Our NLTE 
calculations are always 'restricted', in the sense that the 
temperature and electron density are adopted from the 
LTE structure calculations. 

The examples discussed in this paper show how the 
different ingredients involved in modeling spectral line for- 
mation can begin to be unraveled when high-quality ob- 
servations of spatially resolved solar line profiles are com- 
pared with detailed calculations. High quality solar ob- 
servations with low spatial resolution are still missing at 
most wavelengths and can play a crucial role guiding the- 
ory in the quest for accuracy in the interpretation of stellar 
spectra. 
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